Monitoring of band gap and magnetic state of graphene nanoribbons through 

vacancies 
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Using first-principles plane wave calculations we predict that electronic and magnetic properties 
of graphene nanoribbons can be affected by defect-induced itinerant states. The band gaps of 
armchair nanoribbons can be modified by hydrogen saturated holes. Defects due to periodically 
repeating vacancy or divacancies induce metallization, as well as magnetization in non-magnetic 
semiconducting nanoribbons due to the spin-polarization of local defect states. Antiferromagnetic 
ground state of semiconducting zigzag ribbons can change to ferrimagnetic state upon creation of 
vacancy defects, which reconstruct and interact with edge states. Even more remarkable is that 
all these effects of vacancy defects are found to depend on their geometry and position relative to 
edges. It is shown that these effects can, in fact, be realized without really creating defects 
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I. INTRODUCTION 

Its unusual electronic energy band structure and 
charge carriers resembling massless Dirac Fermions have 
made graphene honeycomb structure an active field of 
researchi^^. Quasi ID graphene ribbons have even 
more interesting electronic and magnetic properties de- 
pending on their size and symmetry 6 . These are edge 
states of zigzag ribbons with opposite spin polarization 7 
and band gaps varying with the width of the ribbon 8 -. 

Theoretical studies have predicted that energetic elec- 
trons and ions can induce polymorphic atomic defects 
, such as vacancies in graphene 9 . Using high-resolution 
TEM the observation of vacancies have been reportecU -. 
Recent studies have shown that vacancies created on two- 
dimensional (2D) graphene by high-energy electron or 
ion irradiation can induce magnetism in a system con- 
sisting of only sp-electron a 11 ! 12 . It has been argued that 
Stoner magnetism with high Tq originates from the spin- 
polarized extended states induced by the vacancy defects, 
while RKKY coupling is suppressed. These effects of de- 
fects on ID semiconducting graphene nanoribbons should 
be more complex and interesting, because their band gap, 
magnetic state and symmetry are expected to intervene. 

Present study based on extensive first-principles as well 
as empirical tight-binding calculations (ETB) has shown 
that the band gap and magnetic state of any armchair or 
zigzag nanoribbons can be modified by single or multi- 
ple vacancies (holes) . The effects of these defects depend 
on their symmetry, repeating periodicity and positions. 
Even if the formation of periodic defects of desired sym- 
metry may not easily be achieved, similar effects can be 
created through the potential difference across a ribbon 
applied by periodically arranged tips. When combined 
with various properties of nanoribbons these results can 
initiate a number of interesting applications. We believe 
that our results are important for further studies, since 
the graphene ribbons can now be produced with preci- 
sion having widths sub 10 nmiS and nanodevices can be 



fabricated thereof^. 



II. MODEL AND METHODOLOGY 

We have performed first-principles plane wave calcu- 
lations within Density Functional Theory (DFT) using 
PAW potentials^. The exchange correlation potential 
has been approximated by Generalized Gradient Approx- 
imation (GGA) using PWQl-i 6 - functional both for spin- 
polarized and spin-unpolarized cases. All structures have 
been treated within supercell geometry using the peri- 
odic boundary conditions. A plane-wave basis set with 
kinetic energy cutoff of 500 eV has been used. In the 
self-consistent potential and total energy calculations the 
Brillouin zone (BZ) is sampled by (1x1x35) special k- 
points for ribbons. This sampling is scaled according to 
the size of superlattices. All atomic positions and lattice 
constants are optimized by using the conjugate gradient 
method where total energy and atomic forces are min- 
imized. The convergence for energy is chosen as 10~ 5 
eV between two steps, and the maximum force allowed 
on each atom is less than 0.02 eV /A. Numerical plane 
wave calculations have been performed by using VASP 
packag o 17 ! 18 



III. VACANCY DEFECTS IN NANORIBBONS 

We start by summarizing the electronic properties of 
graphene nanoribbons which are relevant for the present 
study. We will consider hydrogen terminated nanorib- 
bons if it is not stated otherwise. Armchair graphene 
nanoribbons, AGNR(A) (N being the number of carbon 
atoms in the primitive unit cell) , are non-magnetic semi- 
conductors. The band gap&i 9 ., Eg, of a bare or hydrogen 
terminated AGNR(A^) is small for N = 6m — 2 (m being 
an integer) , but from N = 6m to N = 6m + 2 it increases 
and passing through a maximum it becomes again small 



at the next minimum corresponding to TV = 6m + 4. As 
Ec oscillates with N its value shall decrease eventually to 
zero as N — > oo. Bare and hydrogen terminated zigzag 
nanoribbons, ZGNR(A), are also semiconductors with 
Eg decreasing consistently as N increases for N >8, but 
their edge states give rise to antiferromagnetic (AFM) 
ground state!. We will show that these magnetic and 
electronic properties can be modified by defects originat- 
ing from vacancies or holes created in those ribbons. 

A hole can be created when carbon atoms at the cor- 
ners of any hexagon of an armchair ribbon are removed 
and subsequently remaining six two- fold coordinated car- 
bon atoms are terminated with hydrogen atoms. Fig- 
ure Q] (a) - (b) shows that the electronic structure of 
AGNR(34) is strongly modified by such a hole which is 
placed at the center of the ribbon. The hole repeats itself 
at each supercell which comprises six primitive cells cor- 
responding to a repeat period of 1=6. This ribbon having 
a periodic hole (or defect) is specified as AGNR(iV;Z) and 
has non-magnetic ground state. Despite large separation 
of periodic defect which hinders their direct coupling such 
a strong modification of the band gap is somehow unex- 
pected. However, it is an indirect effect and occurs since 
the itinerant (Bloch) states of band edges are modified 
by the defect. At the end, the direct band gapi^ at the 
T-point has widened from 0.09 eV to 0.40 eV due to a 
defect situated at the center of the ribbon. AGNR (36) 
exhibit the same behavior; namely its band gap increases 
when a similar hole repeats itself at each supercell con- 
sisting of six primitive cell. However, in contrast, the 
band gap of AGNR(38), which is normally larger than 
that of AGNR(34), is reduced if the same hole is intro- 
duced at its center. In addition, states localized around 
the defect have formed flat bands near the edge of valence 
and conduction bands because of their reduced coupling. 

Moreover, the effect of this defect is strongly dependent 
on its position relative to the both edges of the nanorib- 
bon as depicted in Fig. [2] (a). As shown in Fig. [5] (b), 
the changes in band gap depend on N of AGNR (in the 
family specified as 6m + q, q being -2, and +2) as well 
as on the position of the hole. We note that the varia- 
tion of the band gap with the position of the hole relative 
to the edges of the ribbon show similar trends for both 
AGNR (34,6) and AGNR (36;6). However, AGNR (38;6) 
displays an opposite trend. 

ETB calculations indicate that a similar behavior is 
also obtained in AGNR(A;Z)'s having relatively larger 
N (i.e. A=58, 60, 62) in Fig H (c). As expected the 
effect of the hole on the band gap depends on the repeat 
periodicity. As shown in Fig [5] (d), the effects of defect 
decrease with increasing repeat periodicity I. We found 
that larger holes with different geometry and rotation 
symmetry can result in diverse electronic structure and 
confined states, ft should be noted that a repeating hole 
can also modify the mechanical properties. For example, 
the stiffness of a ribbon is reduced by the presence of a 
hole. The force constant, k = d 2 Ex/dc 2 (c being the 
lattice constant) calculated for AGNR(34;6) with a hole 
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FIG. 1: (Color online) (a) Energy band structures of 
AGNR(iV=34;Z=6) with and (b) without a hole consisting 
of six carbon vacancies, (c) Charge density isosurfaces of se- 
lected states. Carbon atoms (represented by black circles), 
which have coordination number lower than 3 are terminated 
by hydrogen atoms (represented by small gray circles). 



at its center (k = 6.03 eV/A) is found to be smaller than 
that without a hole « =7.50 eV/A). 

We now show that different types of vacancies in the 
same armchair nanoribbon gives rise to different changes 
in electronic and magnetic properties. A divacancy cre- 
ated in AGNR(22) can cause to a dramatic change in the 
electronic state of the ribbon when it is repeated with the 
periodicity of 1=5. Such an armchair ribbon is specified 
as AGNR(22;5). The divacancy first relaxes and forms 
an eight fold ring of carbon atoms which is adjacent to 
six hexagons and two pentagons. In Fig. [3] (a) we see 
that the non-magnetic and semiconducting AGNR(22) 
with band gap of Eq=0.18 eV becomes a non-magnetic 
metal, since a flat band derived from the defect occurs 
below the top of the valance band edge and causes to a 
metallic state. 

The effect of a single carbon vacancy becomes even 
more interesting. A single vacancy created in AGNR(22) 
is relaxed and the three-fold rotation symmetry is broken 
due to Jahn- Teller distortion. At the end, a nine-sided 
ring forms adjacent to a pentagon as shown in Fig. [3] 
(b) . Owing to the spin-polarization of the sp 2 -dangling 
bond on the two-fold coordinated carbon atom and ad- 
jacent orbitals at the defect site the system obtains an 
unbalanced spin. In fact, the difference of total charge 
density corresponding to states of different spin states, 
i.e. Apx = Pt ~ Pt 1S non_z ero and exhibits a distribu- 
tion shown in Fig. [3] (b). Because of unbalanced spin, 
AGNR(22;5) gains a net magnetic moment of n=l hb per 



3 




Defect Positions- 



(b) 



00 



(d) 



a 



AliNR (34;6)+6V 
AGMR (36;6)+6V 
I ACiiNR (38;6)+6V 




V AGNR (5»;6)+fiV 
□ ACINR (60;6)+ev 
O AGNR (62;6}+6V 



s. 



T ACJJNR <34:1)+6V 
y AtiPMR (5S:1)+6V 

O A<iNR <62:1)+6V 



12 3 
Defect Position 



(112 3 
Defect Position 



4 8 12 16 
Defect Seperatton "/" 



FIG. 2: (Color online) (a) Positions of a hole in the ribbon 
are indicated by numerals. The nanoribbon without a hole 
is specified by "0". (b) and (c) Variations of band gaps of 
AGNR(iV,Z) with the position of the hole specified as 6V . (d) 
Variation of band gap with repeat periodicity, I. Results for 
N < 58 are obtained by first-principles calculations. 



cell. This attributes a magnetic state to the nanoribbon, 
which was nonmagnetic otherwise. 

Furthermore, the spin degeneracy of some bands re- 
lated with this defect is broken and spin-up and spin- 
down bands are split. The dispersive, non-magnetic band 
at the edge of the valence band becomes partially emp- 
tied, since its electrons are transferred to the fiat spin 
bands below. Eventually, the semiconducting ribbon be- 
comes metallic. 

Not only armchair, but also zigzag nanoribbons are 
strongly affected by defects due to single and multiple va- 
cancies. When coupled with the magnetic edge states of 
the zigzag nanoribbons the vacancy defect brings about 
additional changes. The magnetic state and energy band 
structure of these ribbons depend on the type and ge- 
ometry of the defects. In Fig. 0] (a)-(c), the effect of a 
defect generated from the single vacancy with a repeat 
periodicity of Z=8 is examined in ZGNR(14;8) for three 
different positions. The total energy is 0.53 eV lowered 
when the defect is situated at the edge rather than at the 
center of the ribbon. 

ZGNR(14;8) has a net magnetic moment of /x=1.94 [Ib 
when the defect is situated at the center of the ribbon and 
hence its antiferromagnetic ground state has changed to 
ferrimagnetic state through the magnetic moment of the 
single vacancy. Otherwise, ^ becomes zero, when the 
position of the defect deviates from the center. For ex- 
ample, in Fig. 2] (b) and (c) the sum of the magnetic 
moments of the edge states is not zero, but the net mag- 
netic moment per unit cell becomes zero only after the 
spins of the defect is added. Even if the net magnetic mo- 
ment /x=0, ZGNR(14,8) does not have an antiferromag- 
netic ground state due to the presence of single vacancy. 
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FIG. 3: (Color online) (a) Metallization of the semiconducting 
AGNR(22) by the formation of divacancies with repeat period 
of 1—5. (b) Magnetization of the non-magnetic AGNR(22) 
by a defect due to the single carbon atom vacancy with the 
same repeat periodicity. Isosurfaces around the vacancy cor- 
responding to Apr; namely the difference of the total charge 
density of different spin directions. Solid (blue) and dashed 
(red) lines are for spin-up and spin-down bands; solid (black) 
lines are nonmagnetic bands. 



The edge states, each normally having equal but opposite 
magnetic moments, become ferrimagnetic when a defect 
is introduced. The total magnetic moment of the super- 
cell vanish only after the magnetic moment of defect has 
been taken into account. Since the spin-degeneracy has 
been broken, one can define E G and Eq for majority and 
minority spin states. Not only the magnetic state, but 
also the band gap of zigzag ribbons in Fig. 2] are affected 
by the symmetry and the position of the defect relative to 



4 



(a) 

i.i.t.t.t.t.t.t.i 

*, A.. 

•.•.•.*.•♦*.*.♦. 

E =- 1069.13 eV, JJ.= 1.94 Lt B 
Ej=0.29eV, E*=0.06cV 

(b) 





• I ..... in 



E =-1069.30 eV, H=0|.l B 
eJ=0.21 eV, E ; =0.03 eV 



(C) 

MMMM 




.4 A , i 

E =-1052.82 cV, H=0 ^ B 
E*-O.05cV, Ej=0.O5eV 

( e ) 

. xcv 

4 1 i 

E =-1052.18 eV, |i=3,91|U B 
E^=0.19cV, E*=0.20cV 

(0 

A»4.4j4.4.4 r 4,A 





E = -1069.66 eV. (1=0 (J B 
E?-0.25 eV. E'-O.I5eV 



> • • * - * ■ 

E =-1059.43 cV, M=0|J. B 
E r ,= t>-lt*eV, E j =0.07cV 



FIG. 4: (Color online) Vacancy and divacancy formation in 
an antiferromagnetic semiconductor ZGNR(14) with repeat 
period of 1=8. Calculated total energy, Et (in eV/cell), net 
magnetic moment, fj, (in Bohr magneton /is/cell), band gap 
between spin-up(down) conduction and valence bands, Eq^) 
are shown for each case. Blue and yellow isosurfaces corre- 
sponds to the difference of spin-up and spin-down states 



edges. In Fig. 0] (d)-(f) two defects associated with two 
separated vacancy and a defect associated with a relaxed 
divacancy exhibit similar behaviors. 

An important issue to be addressed here is the break- 
down of Lieb's theorem-- for those zigzag ribbons. Ac- 
cording to Lieb's theorem, the net magnetic moment per 
cell is determined with the difference in the number of 
atoms belonging to different sublattices, and it shall be 
either \i=\ fig or 2 /ib for the cases in Fig. 2] None of the 
cases in Fig. @]is in agreement with Lieb's theorem. Here 
one can consider two features, which may be responsible 
from this discrepancy. First is the strong Jahn- Teller 
distortion and relaxation of carbon atoms at the close 
proximity of the defect. As a result some dangling sp 2 - 
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FIG. 5: (Color online) (a) Monitoring of band gaps Eg by 
applying a local bias voltage Vup across the ribbon, (b) Eg 
versus Vu p applied at the center of AGNR(7V;6) for JV=34, 
36 and 38. (b) E G versus V tip of AGNR(34) for different 
tip positions schematically described at the top. The tip (or 
electrodes) are situated at one of the position 1-3. The repeat 
period is 1=20. 



bonds reconstructed to form new C-C covalent bonds. 
The lowering of the total energy, a driving force for such 
reconstruction, is as high as 0.5-0.6 eV/cell. Second is 
the interaction with the magnetic edge states, which be- 
comes effective for narrow ZGNR's. 

Finally, we note that introducing periodic vacancies or 
defects on a given ribbon appear to be difficult with the 
state of the art technologies. Here we propose a method 
as described in Fig. O to achieve the formation of pe- 
riodic local defects like holes or vacancies. The sharp 
electrodes like STM tips are situated at desired locations, 
such as one of the cases 1-3 in Fig. [5] (a) on the graphcnc 
with a given repeat periodicity. A potential difference, 
Vup common to all electrodes (tips) is applied between 
the tip and underlying insulator through the graphene. 
This way the electronic potential of graphene atoms just 
below the tip is locally lowered or raised depending on 
the polarity of Vu p . Here, the effect of locally and pe- 
riodically applied potential difference has been modelled 
by ETB, where the on-site energies of carbon atoms be- 
low the tip have been changed accordingly. Although 
the present model is crude, it still allows us the realiza- 
tion of monitoring of the properties of nanoribbons. In 
Fig. [5] (b) the variation of Eq with Vn p is calculated 
for AGNR(N;20) with iV=34, 36 and 38 by using period- 
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ically located tips at the center of the ribbon. Because 
of the electron-hole symmetry in the ribbons, the band 
gap variation depends on the magnitude of the bias volt- 
age. In Fig. [5] (c) the variation of the band gap of 
AGNR(42;20) with Vu p and position of the tip is shown. 
In spite of the fact that the modifications of the band 
gaps are not the same as in Fig. [2J the available param- 
eters, such as Vup, I, tip-geometry and its position make 
the monitoring of the properties possible. 



have been found to be identical with the TEM images 
reported earlier 10 . Two different electronic states are 
distinguished: These are: (i) itinerant Bloch states per- 
turbed by defects, (ii) defect induced states. While the 
former is dispersive, the latter give rise to flat bands. The 
reconstruction and spin-polarization of the orbitals at the 
close proximity of the defect give rise to net magnetic mo- 
ments, which, in turn, changes the magnetic ground state 
of the defect-free ribbon. 
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